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Abstract: We have implemented in the RAPGAP program a previously derived sub- 
traction method for next-to-leading order (NLO) corrections in Monte-Carlo (MC) 
event generators, and we show results for jet production in deep inelastic scattering. 
At small Xbj, NLO corrections are comparable to the LO cross section, because of the 
large gluon density. We devise a jet observable that is particularly sensitive to the 
treatment of parton kinematics, which our method is intended to treat correctly. We 
compare the results of the calculation with LO calculations and with the previously 
existing treatment of NLO processes in RAPGAP. Substantial corrections to the event 
distribution are found. 
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1. Introduction 

In contrast to "fixed order" calculations, Monte-Carlo (MC) event generators attempt 
to predict the bulk of the exclusive components of perturbative QCD subprocesses and 
to give the details of the final states. While this is at the expense of model dependence 
in the non-perturbative hadronization, it is also the case that MC event generators 
incorporate a useful approximation to the resummation of large logarithms in their 
parton showering. This provides a successful phenomenology of infrared sensitive cross 
sections. 

Now, it is quite hard to treat non-leading corrections in MC event generators. How- 
ever, this is also quite essential. For example, in deep inelastic scattering (DIS) at small 
Xbj or in diffractive DIS, the gluon density is substantially larger than the quark densi- 
ties, and the next-to-leading order (NLO) gluon-induced process is comparable to the 
leading order (LO) quark-induced process. Therefore, to get a sensible phenomenology, 
the gluon-induced NLO corrections must be included. 

Precisely because MC event generators treat exclusive components of the cross 
section, standard "fixed order" calculations of hard scattering coefficients are not well 
adapted for use in event generators. Some symptoms of this are that in LO, the co- 
efficients contain (5-functions of parton kinematics, and that in NLO the coefficients 
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are singular distributions (generalized functions). They only give results that can be 
compared with hadronic level cross sections after an integration over a wide range of 
parton kinematics, which provides a kind of infrared safety. Standard fixed-order calcu- 
lations approximate the partons as on-shell with zero transverse momentum. For cross 
sections that do not include a wide integration over parton kinematics, only a correct 
treatment of parton kinematics allows a correct conversion of the generalized functions 
to ordinary functions of kinematic variables. In fact, the standard approximations used 
in hard-scattering coefficients cause complications and difficulties in the derivation of 
the algorithms used in event generators. 

In this paper, we will examine the rapidity and transverse momentum distribution 
of the 1-jet inclusive cross section in DIS. Even though this is a jet cross section, it 
is not infra-red safe because its LO hard scattering coefficient is a (5-function, and 
in the approximation that parton virtuality and transverse momentum are ignored, 
the jet momentum can be calculated from the measured 4-momentum of the scattered 
electron. Thus the 5-function appears in the hadron-level jet cross section, in this simple 
parton-model approximation. Similarly, the conventional NLO corrections are singular 
at the parton-model values. Because we only integrate over a narrow range of parton 
kinematics around these singularities, it is not clear how to apply these conventional 
NLO hard scattering coefficients to this cross section. However, it is crucial for a good 
NLO MC method to consistently treat this kind of cross section with NLO accuracy. 

As we argued in [jl], |^, ^, a proper treatment of higher order (NLO and beyond) 
corrections in event generators requires (a) that the parton kinematics should be treated 
more exactly and (b) that the hard scattering coefficients should be ordinary functions, 
not generalized functions. This is in contrast to other proposals @| for matrix element 
corrections to event generators. In [|l| our approach is proposed in general and is 
designed to be generalizable to all orders in perturbation theory. 

More concretely, in this method was applied to DIS in a form specifically de- 
signed for use in event generators using the Bengtsson-Sjostrand algorithm for 
showering. Such event generators include RAPGAP and PYTHIA 0. We have 
now implemented^ this method in RAPGAP, a widely used event generator for both 
inclusive and diffractive DIS. The implementation is based on earlier work by Schilling 

B- 

In this paper we show the results of the 1-jet inclusive cross section calculated 
by our subtraction method. We also compare our method with a method previously 

^Only gluon-induced NLO corrections are implemented, as in The quark-induced NLO terms 
have the compHcation of a soft divergence, which we plan to deal with in a later paper. Gluon-induced 
NLO corrections are important in DIS, because they are not suppressed compared to the LO terms 
whenever the gluon distribution is large. 
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implemented in RAP GAP, which uses a kt-cut apphed to the simple unsubtracted NLO 
cross section. The kt-cut method is intended mainly for the di-jet rate calculation of 
high pt jet events. 

We also compare our results to those of a LO calculation from RAPGAP, where 
the 5-function in the 1-jet inclusive cross section is smeared out only by the parton 
showering and the hadronization. At small Xbj the gluon-induced NLO contributions 
are enhanced, and we expect the cross section with NLO accuracy to be a very broad 
distribution, and the three methods should give very different results. At large Xbj the 
NLO corrections are suppressed by 0{as), we expect to see similar results by these 
methods. 

The organization of this paper is the following. In Sec. ^ we briefly describe the 
subtraction method and the kt-cut method in RAPGAP. In Sec. ^ we show and discuss 
the results of the 1-jet inclusive cross sections for DIS at small Xbj and large Xbj predicted 
by the three different methods. In Sec. |, we discuss the possible applications and future 
work of the NLO RAPGAP with subtraction method. 



2. NLO hard scattering 

In 0, we devised a subtraction method for NLO corrections, and obtained the coef- 
ficient function for the theoretically simple yet phenomenologically significant gluon- 
induced correction to DIS. It was assumed that the Bengtsson-Sjostrand (BS) algo- 
rithm is used for the initial state parton shower, as in RAPGAP. A point-by-point sub- 
traction in the phase space ensures that the hard-scattering coefficient is non-singular 
0. The pdfs needed are specific to the showering algorithm and those needed to 
match the BS algorithm were calculated in 0. The quark-induced NLO corrections 
have the complication of a soft divergence, and we plan to treat this case, matched 
to the BS algorithm, in a later paper; the quark-induced term is a genuine correc- 
tion, order ag compared with the LO process, whereas the gluon-induced term is not 
suppressed at all when the gluon density is large. 

In RAPGAP, the previously implemented method for NLO was the kt-cut method, 
designed mainly to give the correct di-jet rate at large Q^. Here the unsubtracted 
parton-level matrix elements of the NLO processes are used with a cut on the transverse 
momentum kt of the out going partons in the center of mass frame of 7* -|- parton 
scattering, kt cut is required to be large enough so that the NLO contributions are 
smaller than the total cross section. Usually fc^cut = 4GeV^. This gives us a large 
logarithm log{k^ ^^^^/Q"^) in the integrated NLO cross section when is large. 
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3. Rapidity and Transverse Momentum Distribution of DIS 1- 
jet Inclusive Cross Section 



3.1 Theoretical and experimental considerations 

To probe the differences between different treatments of NLO corrections in an event 
generator, we need an infrared sensitive cross section that potentially gives us the max- 
imum difference between the different methods. We choose the rapidity and transverse 
momentum distribution of the 1-jet inclusive cross section, (fa/{dfjdxt) in the labo- 
ratory frame, with the rapidity and transverse momentum measured relative to the 
simplest parton-model values. At the LO parton level, the scattered quark has a trans- 
verse momentum of Ef^ and pseudorapidity r]^^, both of which can be calculated from 
the momentum of the scattered electron, and so we define 

^jet 

If we ignore parton virtuality and transverse momentum, the LO cross section is a 
simple 5-function 

^oc^(..-lMM). (3.3) 

In real QCD, with effects from higher order QCD subprocesses and the hadronization, 
the 5-function will be smeared out. The shape of the smeared 5-function predicted 
by a MC event generator will provide us information about its treatment of the NLO 
corrections. 

We calculate the cross section from the event generator RAPGAP by three meth- 
ods: 

1. LO hard scattering, where jet structure is generated exclusively by parton show- 
ering and hadronization. 

2. With NLO corrections implemented by the kt cut method. 

3. With NLO corrections implemented by our subtraction method. 

Due to the enhanced gluon density at small Xbj, the cross section in the small region 
is likely to give us the maximum difference between the different methods. At large 
Xbj the NLO corrections are suppressed by ct^, and therefore in this region the cross 
section is not sensitive to the difference between the NLO methods. So in this paper 
we calculate this cross section in two different regions, 10-^ < Xbj < 10-^ where we 
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expect to see very different results from tlie tliree calculations, and 0.02 < Xbj < 1, 
where similar results are expected. 

Our calculations are performed under the following conditions, with hadron kine- 
matics appropriate to the collider experiments at HERA: 

• The positron is moving in the —z direction with PI = —27.5 GeV. 

• The proton is moving in the +z direction with = 820 GeV. 

• 20 < g2 < 50 GeV^. 

• In LO and /cj-cut calculations, the LO pdfs Cteq4L are used. 

• In RAPGAP with our subtraction method, we use the BS-algorithm specific pdf 
calculated by a scheme change from the MS pdf, Cteq4M. 

• The jet-finder is PXCONE, with i^^-mode and cone radius= 1. 

• The running scale of the pdfs and as is Q. 

• Initial state showering is turned on, with ordering. 

• Final state showering is on. 

• Hadronization and proton remnant are required. 

• The LO subprocess and (for the NLO calculations) the gluon-induced NLO sub- 
processes are chosen. 

• Only light quarks are considered. 

Figs. D, H and ^ give the results for 10~^ < Xbj < 10~^, while Figs. ^ and ^ give the 
results for 0.02 < x^j < 1. 

3.2 Cross Section at small Xbj 

In Fig. 1^ we plot d?a / {dfidxt) as a function of Xt and fj for the small Xbj region, 10""^ < 
Xbj < 10~^. It clearly shows that the LO (graph c) prediction is strongly peaked at the 
PM value Xt = 1 and = 0, as expected. The cross section predicted by the fct-cut 
method (graph b) is somewhat broader but still peaked at the PM values. In contrast, 
the prediction from our subtraction method (graph a) is broadly distributed; although 
it has a maximum close to the PM position, the peak is not at all 5-function like. This 
indicates that when the gluonic NLO term is comparable to the LO term, a correct 
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(a) RAPGAP-sub (b) RAPGAP-k^-cut (c) RAPGAP-LO 




'^X 



Figure 1: Transverse momentum and rapidity distribution of 1-jet inclusive cross section 
for 10^'' < Xbj < 10~2. (a) RAPGAP with the subtraction method, (b) RAPGAP with the 
kf-cut method and (c) LO RAPGAP. Note that the vertical scales for these three graphs are 
different. 

treatment of off-shell parton kinematics has a substantial effect on the shape of the 
cross section. 

For a more quantitative comparison between these methods, we now concentrate 
on the one parameter cross sections da/dr) and da/dxt, shown in Figs. |^ and ^. In 
each figure, we show in graphs (a), (b) and (c) the components of the 1-jet inclusive 
cross section separated according to the number of jets in the final state: 1 jet, 2 jets 
and > 3 jets. Then in graph (d), we give the total 1-jet inclusive cross section. The 
solid curves, dashed curves and the dotted curves are the predictions of our subtraction 
method, the kt-cut method, and the LO calculation, respectively. 

The differences between the predictions are evident. The subtraction method gives 
the broadest distribution, while the kt-cut method and the LO results show progres- 
sively more strongly peaked distributions. A comparison with data should fairly easily 
test our claim that the subtraction method is more correct. We also see that the num- 
ber of 1-jet events and the number of 2-jets events predicted by our subtraction method 
are of the same order, while the fct-cut method and LO predict more dominance of 1-jet 
events^. Also note that, in Fig. ^ the peaks of the curves are not at the PM position 
fj = 0. Even though the peak corresponds to a smeared out ^(17), the peak has shifted 
to f] around 1/2. The reason is that showering and hadronization of the jet changes a 

^The number of events is roughly proportional to the area under the curve divided by the number 
of jets in an event. 
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Figure 2: Rapidity distribution of 1-jet inclusive cross section for 10~^ < Xbj < 10~^. The 
four graphs show the following contributions: (a) 1-jet events, (b) 2-jet events, (c) > 3-jet 
events, (d) total cross section. The solid curve, dashed curve and dotted curve are predictions 
from the subtraction method, the kt-cut method and the LO RAPGAP, respectively. 



massless on-shell quark into a massive quark, and the kinematics of 2-body collisions 
shows that this moves the rapidity of the quark towards the proton's rapidity. Since 
the proton is chosen to move in the +z direction, this means that the rapidity of the 
jet is biased to a positive value relative to the simple parton model value. 

3.3 Cross Section at large Xbj 

In Figs. ^ and ^ we show the results for large Xbj. The parameters are the same as 
before, except that we choose to analyze events with 0.02 < Xbj < 1. In contrast to the 
small Xbj results, the three different methods give quite similar results, all peaked near 
fj = and Xt = 1. In addition, most of the events are 1-jet events. This is because 
gluon distribution is now much smaller, and the LO process dominates. 
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Figure 3: Transverse momentum distribution of 1-jet inclusive cross section. The coding of 
the curves is the same as in Fig. |2[ 



4. Applications and Future Work 

We have implemented our subtraction method in RAPGAP for DIS processes. The 
subtraction method is designed to give perturbative QCD predictions for the hard scat- 
tering with a correct treatment of parton kinematics. We chose to calculate an infra-red 
sensitive observable (Pa/dxtdfj using three different methods, and to perform the calcu- 
lation in a region where the NLO correction is large because of a large gluon density. We 
see large differences between the subtraction method and the kt-cvX method, thereby 
demonstrating the importance of the treatment of parton kinematics. 

There is no experimental analysis on (Pa/dxtdfj available yet, current analyses 
lOl , [Till being in the Breit frame. Since our subtraction method, combined with the 



parton shower, is intended to give more accurate perturbative QCD predictions for 
infrared sensitive cross sections, it would be interesting to do the experimental analysis 
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Figure 4: Rapidity distribution of 1-jet inclusive cross section for 0.02 < Xbj < 1- The 
coding of the curves is the same as in Fig. 0. 



of the jet cross section in the laboratory frame relative to the position of parton-model 
jet, and to compare it with the MC calculations using our subtraction method. 

In the future we would like to work on two directions. One is to implement the 
subtraction method for difFractive DIS where there is also a large gluon distribution. 
Another is to work on the quark-induced NLO subprocess in order to get a complete 
NLO MC event generator. 
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Figure 5: Transverse momentum distribution of 1-jet inclusive cross section for 0.02 < Xbj < 
1. The coding of the curves is the same as in Fig. |^. 
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